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1. Introduction 

The confinement mechanism in Coulomb gauge QCD has received a lot of attention recently. 
It is firstly discussed by Gribov in '70s that the instantaneous interaction provides the long range 
interaction [jlj], and this is further elaborated by Zwanziger recently Zwanziger showed that 
the color-Coulomb potential which is the instantaneous interaction energy between heavy quarks 
is stronger than a physical potential. This inequality tells us that the necessary condition for the 
physical potential being a confining potential is that the color-Coulomb potential is also a confining 
potential, i.e., "no confinement without color-Coulomb confinement" ph. 

The recent Monte Carlo simulations in the SU(2) and SU(3) lattice gauge theories showed 
that the color-Coulomb potential rises linearly with distance, and its string tension has 2^3 times 
larger value than that of the Wilson potential, which is an expected result from the Zwanziger's 
inequality. [Q, §J. In addition, it was shown that an asymptotic scaling violation for the color- 
Coulomb string tension is weaker than that of the Wilson string tension [||, ^]. Furthermore, it has 
been shown that the color-Coulomb potential is a confining potential even in the deconfinement 
phase [Q, [50. Thus the color-Coulomb string tension does not serve as an order parameter for the 
confinement/deconfinement phase transition. 

In Gribov-Zwanziger confinement scenario, the near-zero modes of the Faddeev-Popov (FP) 
operator play a crucial role to produce the singular behavior of the ghost propagator in the infrared 
region. The color-Coulomb potential in the color-singlet channel is given by 

v^*-)0 = s 2 Tr[rW|^ (1.1) 

where & is the Green's function of the FP operator and (•) denotes an Euclidean expectation value. 
T" (a = 1, 8) are the generators of su(3) Lie algebra. The singular behavior of the ghost prop- 
agator in the infrared region leads to the long-range interaction of the color-Coulomb potential 
which is responsible for the color confinement. 

Recently Greensite, Olejnfk and Zwanziger discussed the self-energy of an isolated quark and 
derived the necessary condition for the color confinement The authors studied the spectrum of 
the FP operator in Coulomb gauge using SU(2) lattice gauge simulation and confirmed that the nec- 
essary condition is satisfied. In this study, we investigate the distribution of the FP eigenvalues in 
SU(3) lattice gauge simulations, and check whether the necessary condition for color confinement 
is satisfied or not. 

2. Color-Coulomb self-energy 

The Coulomb gauge Hamiltonian can be expressed as the sum of the transverse part and the 
instantaneous part [^] : 

H = l -jd\{{Ef r (x,t)f+B^x,t) 2 ) + l -jd^y J d 3 zp a (y,t)r ah (y,z;A tr )p h (lt)- (2.1) 

Here Ef' tr wee the transverse components of the color electric field, Bf the color magnetic field, 
p a (x,t) the color charge density. The kernel of the instantaneous interaction is given by 

r ab (y,z;A' r ) = J d 3 x^ ac (y,x]A' r )(-Vi)W cb (x,z]A tr ), (2.2) 



2 



Infrared divergence of the color-Coulomb self-energy in Coulomb gauge QCD 



Y. Nakagawa 



where A" ,tr are the transverse components of the gluon field. 5f is the Green's function of the FP 
operator M ab = -dpf = -8 ab d 2 - gf abc A c / r d h On a lattice, the FP operator is an 8V 3 x 8V3 
sparse matrix (V3 is the lattice 3-volume) and expressed in terms of SU(3) spatial link variables £/, 

as 



M$ = £>eTr [{T a ,T b } [Ufa) + Vi(x-t)) 8 x , y 

i 

- 22* W,£*)£^ - 2T a T b Ui{x - i)S yA 



(2.3) 



E c = g 2 C D l dX ^^l , (2.6) 



The color-Coulomb self-energy for an isolated color charge, whose energy diverges in the 
infrared limit in a confining theory, is []8|] 

E c = Tr[T a T b }g 2 (r ab (x,x;A tr )). (2.4) 

The color-Coulomb self-energy is ultraviolet divergent in the continuum limit both in an abelian 
and a non-abelian gauge theories, and can be regularized by introducing the cutoff. The interesting 
point is that the infrared divergence may exist in a confining theory at infinite volume. 
We define the normalized density of the FP eigenvalues 

jV(A,A+AA) 

P(A)= 8V3AA ' (2 " 5) 

where N(X,X + AA) is the number of eigenvalues in the range [A, A + AA]. By expanding the 
Green's function of the FP operator in terms of the eigenvectors <j>" (x) and the eigenvalues A ;1 of 
the FP operator, we find 

rlmax (p(A)F(A)) 

"A2 

where Cd(> 0) is the Casimir invariant for the representation D and the upper limit of the inte- 
gration Amax corresponds to the UV lattice cutoff. In the Gribov-Zwanziger scenario the gauge 
configurations are restricted to the Gribov region, and therefore the lower limit of the integration is 
zero. F n are the expectation values of the negative Laplacian in the FP eigenmodes, 

F n = J d 3 x(t>: a (x)(-V 2 )r n (x)- (2.7) 

If the condition 

hm mpm >0 (2 .8) 

A^O A 

is satisfied in the infinite volume limit, the color-Coulomb self-energy diverges in the infrared 
region. This is the necessary condition for the color confinement [JSj] . 

The FP eigenvalue density of the near-zero modes is closely related to the infrared behavior 
of the color-Coulomb potential. From Eqs. ( |1.1[ ) and (2.4), the color-Coulomb self-energy can be 
expressed as 

d3 P f>/=ft_ f Ad \P\^ar, 







Vc( *-*>= / (2#w)=y i^™ 1 ^- (2 - 9) 



Here we introduce the ultraviolet cutoff A. If the condition (2.8) is satisfied, the color-Coulomb 



self-energy diverges in the infrared region. Accordingly, the right-hand side of (2.9) diverges in 
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the infrared limit. It means that the color-Coulomb potential is more singular in the infrared region 
than the Coulomb potential V(p) ~ 1/|/?| 2 - 

The FP operator is the negative Laplacian in the abelian gauge theory. Thus the FP eigenfunc- 
tions are the plane waves and A = k 2 . By counting the number of states in momentum space, it is 
easy to show that 

p(A) = ^, F(A)=A, (2.10) 



in the infinite volume limit. Obviously the necessary condition (2.8) is not satisfied in this case. 



3. Numerical simulations 

We calculate the FP eigenvalue density by the SU (3) lattice gauge simulations in quenched 
approximation. The lattice configurations are generated by the heat-bath Monte Carlo technique 
with the Wilson plaquette action. In these simulations we adopt the iterative method to fix a gauge. 
We used the ARPACK package to evaluate the lowest 1000 eigenvalues and corresponding eigen- 
vectors of the FP operator. 

3.1 (p(A)F(A))/A in the confined phase 




X X 

Figure 1: (a) The FP eigenvalue density in the confinement phase, (b) The average Laplacian F(X) in the 
confinement phase. 



Figures [l](a) and (b) show p (A) and F(A) at j8 = 6.0 on a variety of lattice sizes. We see the 
accumulation of the near zero modes of the FP operator at larger lattice volume. On the other hand, 
p(A) is almost saturated above A ~ 0.15. F(A) becomes flat at smaller value of A as the lattice 
volume increases and it seems that as A — > the average Laplacian approaches positive constant in 
the infinite volume limit. 

To compare the behavior of p(A) in the non-abelian theory with that in the abelian theory, 
we plot p(A)/vX in Fig. ^(a). In the abelian theory it is constant because p(A) ~ \Pk (see Eq. 
( |2.10 )). In the non-abelian theory, we observe that p(A)/\/A is almost constant above A ~ 1 
[GeV]. In contrast, at small A, p(A)/\/A is not constant and it seems to diverge as A — > in the 
infinite volume limit. The FP eigenvalue density in the non-abelian theory shows a completely 
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different behavior at small A compared to that of the abelian theory and we see the enhancement of 
the near-zero modes of the FP operator. In the Gribov-Zwanziger scenario, these near-zero modes 
cause the color-Coulomb potential to be more singular in the infrared region than the simple pole. 




.4 



Figure 2: (a) p(A)/VA in the confinement phase, (b) p(X)F(X)/X vs. A in the confinement phase. The 
solid curve represents p{X)F(X)/X in the case of the abelian theory, p(X)F{X)/X = \/~X/An 2 . 



In Fig. |2](b) we plot p(A)F(A)/A as a function of A. As A approaches to 0, p(A)7 ? (A)/A 
decreases for the free field (see Eq. ( |2.10| )) while increases for the interacting field. From this 
figure, we expect that p(A)F(A)/A diverges or goes to positive constant, and it is unlikely that 
p(A)F(A)/A goes to zero as A — > in the infinite volume limit. Therefore, we conclude that the 
color-Coulomb self-energy of an isolated color charge is infrared divergent in SU(3) lattice gauge 
theory. 

3.2 (p (A)7 r (A))/A in the deconfinement phase 

p(A), F(X) and p(A)F(A)/A at T /T c ~ 1.5 (T c is the critical temperature of the phase transi- 
tion) are displayed in Figs. ||(a), (b) and (c). These figures show that there are no drastic changes of 
these behaviors in the deconfinement phase. This is consistent with the fact that the color-Coulomb 
potential is not screened above the critical temperature. We conclude that the necessary condition 
for color confinement is satisfied even in the deconfinement phase. 

4. Conclusions 

We have calculated the eigenvalue distribution of the FP operator in Coulomb gauge using 
quenched SU(3) lattice gauge simulations. In the confinement phase, we observe the accumula- 
tion of the near-zero eigenvalues of the FP operator at large lattice volumes. We conclude that 
the confinement criterion is satisfied in the SU(3) gauge theory. Accordingly, the color-Coulomb 
potential becomes more singular than the simple pole in the infrared region. This supports the 
Gribov-Zwanziger confinement scenario. The results we obtained are qualitatively consistent with 
those of the SU(2) lattice simulation carried out by Greensite et al. 

The near-zero modes of the FP operator survive above the critical temperature, and the be- 
haviors of the FP eigenvalue density and the average Laplacian in the deconfinement phase are 
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Figure 3: (a) The FP eigenvalue density in the deconfinement phase, (b) The average Laplacian F(X) in the 
deconfined phase, (c) p (X)F (X) / X vs. X in the deconfinement phase. 



qualitatively the same as in the confinement phase. Accordingly, the confinement criterion is sat- 
isfied even in the deconfinement phase in SU(3) gauge theory. This would indicate that confining 
features survive even in the deconfinement phase, and we expect that further studies in Coulomb 
gauge provide insight into the understanding the strongly correlated quark-gluon plasma. 
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